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Abstract 

We use the deformed sine-Gordon models recently presented by Bazeia et al pQ to 
discuss possible definitions of quasi-integrability. We present one such definition and use 
it to calculate an infinite number of quasi-conserved quantities through a modification 
of the usual techniques of integrable field theories. Performing an expansion around 
the sine-Gordon theory we are able to evaluate the charges and the anomalies of their 
conservation laws in a perturbative power series in a small parameter which describes 
the "closeness" to the integrable sine-Gordon model. Our results indicate that in the 
case of the two-soliton scattering the charges are conserved asymptotically, i.e. their 
values are the same in the distant past and future, when the solitons are well separated. 

We back up our results with numerical simulations which also demonstrate the exis- 
tence of long lived breather-like and wobble-like states in these models. 



1 Introduction 



Solitons and integrable field theories play a central role in the study of many non-linear phe- 
nomena. Indeed, it is perhaps correct to say that many non-perturbative and exact methods 
known in field theories are in one way or the other related to solitons. The reason for that is 
twofold. On one hand, the appearance of solitons in a given theory is often related to a high 
degree of symmetries and so to the existence of a large number of conservation laws. On the 
other hand, in a large class of theories the solitons possess a striking property. They become 
weakly coupled when the interaction among the fundamental particles of the theory is strong, 
and vice-versa. Therefore, the solitons are the natural candidates to describe the relevant 
normal modes in the strong coupling (non-perturbative) regime of the theory. Such relation 
between the strong and weak coupling regimes have been observed in some (1+1) dimensional 
field theories, as, for example, in the equivalence of the sine-Gordon and Thirring models [2J, 
as well as in four dimensional supersymmetric gauge theories where monopoles (solitons) and 
fundamentals gauge particles exchange roles in the so-called duality transformations [3]. 

The exact methods to study solitons in (1 + l)-dimensional field theories involve many 
algebraic and geometrical concepts, but the most important ingredient is the so-called zero 
curvature condition or the Lax-Zakharov-Shabat equation |4J. All theories known to possess 
exact soliton solutions admit a representation of their equations of motion as a zero curvature 
condition for a connection living in an infinite dimensional Lie (Kac-Moody) algebra [5j [6] . 
In fact, in (1 + 1) dimensions such zero curvature condition is a conservation law, and the 
conserved quantities are given by the eigenvalues of the holonomy of the flat connection 
calculated on a spatial (fixed time) curve. On the other hand, many techniques, like the 
dressing transformation method, for the construction of exact solutions are based on the 
zero curvature condition. In dimensions higher than two the soliton theory is not so well 
developed, even though many exact results are known for some 2-dimensional theories such 
as the CP N models [TJ, as well as in four dimensional gauge theories where instanton and 
self-dual monopoles are the best examples [8j. Some approaches have been proposed for the 
study of integrable field theories in higher dimensions based on generalizations of the two 
dimensional methods like the tetrahedron equations [9] and of the concept of zero curvature 
involving connections in loop spaces [TO] . 

Another important aspect of integrable field theories is that they serve as good appro- 
ximations to many physical phenomena. In fact, there is a vast literature exploring many 
aspects and applications of perturbations around integrable models. In this paper we want 
to put forward a technique that, so far as we know, has not been explored yet and which 
suggests that some non-integrable theories often possess many important properties of fully 
integrable ones. We put forward and develop the concept of quasi-integrability for theories 
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Figure 1: Plots of solitons for various values on n 



that do not admit a representation of their equations of motion in terms of the Lax-Zakharov- 
Shabat equation, but which can, nevertherless, be associated with an almost flat connection 
in an infinite dimensional Lie algebra. In other words, we have an anomalous zero curvature 
condition that leads to an infinite number of quasi-conservation (almost conservation) laws. 

Moreover, in practice, in physical situations, like the scattering of solitons, these charges 
are effectively conserved. The striking property we have discovered is that as the scattering 
process takes place the charges do vary in time. However, after the solitons have separated 
from each other the charges return to the values they had prior to the scattering. Effectively 
what we have is the asymptotic conservation of an infinite number of non-trivial charges. 
There are still several aspects of this observation that have to be better understood but we 
believe that if our results are indeed robust then such asymptotic charges could play a role in 
many important properties of the theory like the factorization of the S-matrix. 

We introduce our concept of quasi-integrability through a concrete example involving a 
real scalar field theory in (1 + 1) dimensions which is a special deformation of the sine-Gordon 
model. The scalar field tp of our theory is subjected to the potential 



where n is a real parameter which in the case n = 2 reduces the potential to that of the 
sine-Gordon model, i.e. V (p, 2) = ^ [1 — cos (4 </?)]. 

This potential ( II. ip is a slight modification of that introduced by D. Bazeia et al, [TJ, in 
the sense that we take the absolute value of sin (p to allow n to take real and not only integer 
values. 



V (<p, n) = — tan 2 ip [1 



sin (p 
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Figure 2: Plots of the potential against ip, for three values of the parameter e, 

where n = 2 + e, namely e = —0.5 (fine line), e = 0.0 (thick line), and e = 0.5 (dashed 
line) . 



The potential (11. ip has an infinite number of degenerate vacua that allow the existence 
of solutions with non-trivial topological charges. It is worth noticing that the positions of 
the vacua are independent of n and so they are the same as in the sine-Gordon model, i.e. 
V'vac. = m f , with m being any integer. 

The model with the potential (II .ip is fully topological {i.e. it satisfies its Bogomolnyi 
bound for any n) and so its one soliton field configurations are known in an explicit form. 
They are given by: 



ip = arcsm 



e 2r 



1 + e 



21' 



1/n 



r=± (x -; i 7°>. (i.2) 

v l — f 



where the velocity v is given in units of the speed of light, and the signs correspond to the 
kink (+), and anti-kink (— ), with topological charges +1, and —1 respectively. In Figured] 
we plot the fields of a one soliton configuration for various values on n. We see from this plot 
that the n = 2 case does not appear to be very special; all soliton fields look very similar and 
the solitons for different values of n differ only in their slopes. 

We are going to use this model to back up our discussion of quasi-integrability and so next 
we look at n = 2 + e, with e small. In Figure|2]we plot the potential (II. ip for e = —0.5; 0.0; 0.5. 
Of course, in this case the kinks solutions are given by (II. 2p with n replaced by 2 + e. In 
Figure [3] we plot the one kink solutions (ll.2j) for the potentials shown in Figure [21 i. e. for 
e = —0.5; 0.0; 0.5. Note that they connect the vacua ip vac = to ip vac , — |, as x goes from 
— oo to +oo, with the slope of the kink increasing as the value of e decreases. 

In this paper we study the concept of quasi-integrability in the context of the theory (11.11) 
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Figure 3: Plots of the kink solutions (| 1 . 2 j) against x, with t = and xq = 0, for three 
values of the parameter e, where n = 2 + e, namely e = —0.5 (fine line), e = 0.0 (thick 
line), and e = 0.5 (dashed line). 



from the analytical and numerical points of view^. Our approach and the main results of this 
paper can be summarised as follows (more details are given in the following sections): 

We first consider a real scalar field theory with a very general potential V (</?), and construct 
a connection based on the si (2) loop algebra which, as a consequence of its equations of 
motion, satisfies an anomalous zero curvature condition. Using a modification of the methods 
employed in integrable field theories we construct and infinite number of quasi-conserved 
charges for such a theory, i.e. 

= ~ « (2n+1) (t) n = 0, ±1, ±2, . . . , (1.3) 

where the anomalies cr 2n+1 ) (t) are non-zero due to the non-flatness of the connection A^. 
The charges Q^ 1 ^ are in fact conserved, i.e. a^ 1 ^ = 0, and linear combinations of them 
correspond to the energy and momentum. 

We then restrict ourselves to the case of the potential ( 11.11) and set up a perturbative 
expansion around the sine-Gordon theory. We expand all the quantities, equations of motion, 
field tp, charges and anomalies, in powers of the parameter e, related to n appearing in (11. ip 
by n = 2 + e. For instance, we have 

Q (2n+1) = g (2n+l) + £ Q (2n + 1) + Q ^ ? a ( 2n+1) = £ + Q ^ (L4) 

The anomalies vanish in the lowest order (order zero) in e because they correspond to the 
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sine-Gordon theory which is integrable and so all the charges Q n are conserved during the 
dynamics of all field configurations. 

In this paper we concentrate our attention on the evaluation of the first non-trivial charge 
and its anomaly, namely and a^ 3 \ but our calculations can easily be extended to the other 
charges. We considered the case of the scattering of two kinks and also of a kink/anti-kink in 
the theory (11.11) . where the solitons are far apart in the distant past and future, and collide 
when t ~ 0. We found that the first order anomaly a\ , vanishes when integrated over the 
whole time axis. Therefore, from (I1.3P we see that 

Qf) (t = +oo) = Q? ] (t = -oo) . (1.5) 

Consequently, the scattering of the solitons happen in a way that, to first order in e at least, 
the charge is asymptotically conserved. That is a very important result, and if one can extend 
it to higher orders and higher charges one would prove that effectively the scattering of solitons 
in the theory (11.11) takes place in the sdame way as if the theory were a truly integrable theory. 
We have also analyzed the first order charge for the breather solution of the theory (II. ip . 
and found that even though the charge is not conserved, it oscillates around a fixed value. In 
other words, the first order anomaly vanishes when integrated over a period ^, and so from 
(OD we find that 

Qi ) (t) = Q? (t+l), (1-6) 

where v is the angular frequency of the breather. That means that the period of the charge 
is half of that of the breather. 

We have performed many numerical simulations of the full theory (II -ip using a fourth 
order Runge-Kuta method, and using various lattice grids to make sure that the results are 
not contaminated by numerical artifacts. We have found reliable results with lattice grids of 
at least 3001 points, where the kinks were of size ~ 5 points. Most of the simulations in this 
paper were performed with lattices of 10001 points (i.e. well within this reliability). The main 
results we have found are the following: We have found that if | e | does not get very close 
to unity the kinks and the kink/anti-kink scatter without destroying themselves and preserve 
their original shapes, given in (II. 2p . For small values of | e | the anomaly integrates to zero 
for large values of the time interval, and so the charge is asymptotically conserved within 
our numerical errors. This is an important confirmation of our analytical result described 
above, and is valid for the full charge and not only for its first order approximation as in (II. 5p . 

One of the important discoveries of our numerical simulations is that the theory (II. ip also 
possesses very long lived breather solutions for e ^ 0, which correspond to non-integrable 
models. These long-lived breathers were obtained by starting the simulations with a field 
configuration corresponding to a kink and an anti-kink. As they get close to each other 



5 



they interact and readjust their profiles and some radiation is emitted in this process. We 
absorbed this radiation at the boundaries of the grid and the system stabilized to a breather- 
like configuration. For n = 2 the resultant field configuration was the exact (and analytically 
known) breather while for e small those breather-like fields lived for millions of units of time. 
As one changed e and made it come close to unity the quasi-breathers radiated more and for 
even larger values they eventually died. We also looked at the anomalies for such breather-like 
configurations and have found a good agreement with our analytical results described above. 
The anomaly, integrated in time, does oscillate and for small values of e the charge is periodic 
in time. That is, again, in agreement with the analytical result (II. 6p . Notice however, that 
the numerical result is stronger in the sense that it corresponds to the full charge and not only 
to its first order approximation as in (11.61) . 

We have also performed similar numerical simulations of wobbles [12J which correspond to 
configurations of a breather and a kink. Again, such configurations were obtained by starting 
the simulation with two kinks and an anti-kink. As the three solitons interact and adjust their 
profiles they radiate energy and this radiation had been absorbed at the boundaries of the 
grid. Eventually the system has evolved to a breather and a kink and for small values of e the 
resultant configuration was quite stable, and for n = 2 it agreed with the analytically known 
configuration of a wobble. Again, we believe that this is a very interesting result which shows 
that non-integrable theories can support such kinds of solutions. 

Our results open up the way to investigate large classes of models which are not really 
exactly integrable but which possess properties which are very similar to those of integrable 
field theories. We believe that they will have applications in many non-linear phenomena of 
physical interest. 

The paper is organized as follow: in section [2] we introduce the quasi-zero curvature con- 
dition, based on the sl(2) loop algebra, for a real scalar field theory subject to a generic 
potential, and construct an infinite number of quasi-conserved quantities. In section [3] we 
perform the expansion of the theory (11.11) around the sine-Gordon model, and evaluate the 
first non-trivial charge and its corresponding anomaly. The numerical simulations, involving 
the two solitons scattering, breathers and a wobble, are presented in section HI In section 
[5] we present our conclusions; the details of the sl(2) loop algebra, charge calculations and 
e-expansion are presented in the appendices. 
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2 The quasi zero curvature condition 



We shall consider Lorentz invariant field theories in (1 + l)-dimensions with a real scalar field 
if and equation of motion given by 

dV(ip) 



d 2 ip + 



dip 



0, 



(2.1) 



where V (</?) is the scalar potential. Thus we want to study the integrability properties of such 
theory using the techniques of integrable field theories jH [6] . We then start by trying to 
set up a zero curvature representation of the equations of motion (12. ip . and so we introduce 
the Lax potentials as 



A A 



A_ 



(lo 2 V 



m 



Oi — % to — — Fx 

dip 



■ to d^ip F 



(2.2) 



2 2 

Our Lax potentials live on the so-called sl(2) loop algebra with generators &2n+i and F n , with 
n integer; their commutation relations are given in Appendix [A] The parameters to and m 
are constants, and they play a special role in our analysis. Note, that the dynamics governed 
by (12.11) does not depend upon them, but since they appear in (I2.2p they will play a role 
in the quasi-conserved quantities that we will construct through the Lax equations. In the 
expression above we have used light cone coordinates x± — \ (t ± x), where d± = d t ± d x , and 
d+<9_ = d? -81 = d 2 . 



The curvature of the connection (I2.2p is given by 
d + A--d-A++[A 



F+_ 



with 



X 



% LO 



d-ip 



A_ 



d 2 V 

dip 2 



2 



dV 

dip 



F 



+ lo V — m 



As in the case of the sine-Gordon model where the potential is given by 

^SG = ^ [1 - cos (4^)] 
lo 



(2.3) 



(2.4) 



(2.5) 



we find that X, given by (12 .4p . vanishes when we take lo = 4 and m = 1. Then the curvature 
(12. 3p vanishes when the equations of motion (12.11) hold. The vanishing of the curvature allows 
us to use several powerful techniques to construct conserved charges and exact solutions. We 
want to analyze what can be said about the conservation laws for potentials when X does not 
vanish but can be considered small. 
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In general, the conserved charges can be constructed using the fact that the path ordered 
integral of the connection along a curve T, namely P exp J r da ^f— , is path independent 
when the connection is flat (5J El [13]. Here, we will use a more refined version of this technique 
and try to gauge transform the connection into the abelian subalgebra generated by the &2n+i- 
We follow the usual procedures of integrable field theories discussed for instance in [HI [151 IB] • 
An important ingredient of the method is that our sl(2) loop algebra Q is graded, with n being 
the grades determined by the grading operator d = T 3 + 2 (see appendix [A] for details) 

Q = S ^ J Qn\ [Gm , Qn] C Gm+n ] [d,Q n ] = nQ n . (2.6) 

n 

We perform a gauge transformation 

A^a^gA^g- 1 -d^gg- 1 (2.7) 

with the group element g being an exponentiation of generators lying in the positive grade 
subspace generated by the F n 's, i.e., 



g = exp 



5^ Cn Fn 



Ln=l 



(2.i 



with ( n being parameters to be determined as we will explain below. Under (12.71) the curvature 
(I2.3P is transformed as 

F + _ -> gF + _ g- 1 = d + a_ - d_a+ + [a+ , a_ ] = X g F l g~\ (2.9) 

where we have used the equations of motion (12.11) to drop the term proportional to F in 
(12. 3p . The component A^ of the connection ( 12. 2p has terms with grade and —1. Therefore, 
under ( 12. 7p it is transformed into a_ which has terms with grades ranging from —1 to +oo. 
Decomposing a_ into grades we get from ( 12. 7ft and (12. 8ft that 

a- = X -b„ x (2.10) 

- ~G [b- 1: F 1 }- l -ud^F 

- \t2 [b^,F 2 ] + ^C! [[b- 1 ,F 1 ],F 1 }- l -ud-<pCi [F 1 ,F ]-d^ 1 F 1 



- -Cn[b-i,F n } + .... 

Next we note that one can choose the parameters Cn recursively by requiring that the 
component in the direction on cancels out in a_. Thus we can put that £i = -^ujd-tp, 
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and so on. In the appendix |B] we give the first few C„'s obtained that way. In consequence, 
the component a_ is rotated into the abelian subalgebra generated by the &2n+i- Note that 
this procedure has not used the equations of motion (12.11) . We then have 



6-1 + Yl a 



(2n+l) , 

°2n+l 



n=0 



and the first three components are 



(2.1i; 



,(!) 



,(3) 



,( 5 ) 



(2.12) 



7 



32 



lb 



With the (nS determined this way we perform the transformation of the A + component of 
the connection ( 12. 2p . Since the £ n 's are polynomials of ^--derivatives of ip (see appendix IB"]) 
and since there will be terms involving x + -derivatives of £ n 's, we use the equations of motion 
to eliminate terms involving <9 + <9_</?. Due to the nonvanishing of the anomaly term X in (12. 3p 
we are not able to transform a + into the abelian subalgebra generated by the &2n+i- We find 
that a+ is of the form 



a + 



n=0 



(2n+l) 



6 2 n+l + C 
n=2 



(») 



where 



,(3) 



,(5) 



1 r 2 V 

- w V - m 

2 L 

-a; o_v9- -iujd-(pX, 

4 dip 2 



(2.13) 



(2.14) 



— zu; 3 (<9_</?) 3 X + — u 4 d 2 _<p(d_<p) 2 iud^d 2 _X + -iudl ipd^X 

8 16 ay? 2 2 

2 4 ay? 

with X given in (I2.4p . and V being the potential (see (12. ip ). See appendix |B] for more details, 
including the terms involving c+ . 

The next step is to decompose the curvature ( 12.91) into the component lying in the abelian 
subalgebra generated by &2n+i and one lying in the subspace generated by F n . Since the 
equation of motion (12.11) has been imposed, it turns out that the terms proportional to the 
F n 's in the combination — <9_a+ + [ a + , a_ ] — X g F± g~ l , exactly cancel out. We are then left 



9 



with terms in the direction of the &2n+i only. Therefore, the transformed curvature (12.91) leads 
to equations of the form 

d+a (2n+1) - d^ n+l) = ^ 2n+ ^ n = 0, 1, 2, . . . (2.15) 

with ^" 2,n+l > being linear in the anomaly X given in ( I2.4p . and the first three of them being 
given by 

= 0, (2.16) 

0(3) = iud 2 _(pX, 

1 X. 



iuj 



^u 2 (d-<p) 2 d 2 _ip + di<p 



(2n+l) 



Working with the x and t variables we have that (" 12. 15[) takes the form d t a^ 2n+1 ^ — d x 

<Uyl " +U V" +1 > + af n+1) \T = ^ (2.17) 



dt 2 
with 

/OO POO 
dxaf n+l \ a ( 2 " +1 )=/ dx^ 2n+l \ (2.18) 

-oo J — oo 

As we are interested in finite energy solutions of the theory ( 12. ip we are concerned with 
field configurations satisfying the boundary conditions 

d^ip — > ; V (ip) — > global minimum as x — > ±oo. (2-19) 

Therefore from ( 12. 2 p we see that 

A + ->■ - (w 2 Kac. - 6i, A_ ->■ - 6_j as x ->■ ±oo, (2.20) 

where V vac , is the value of the potential at the global minimum which, in general, is taken to 
be zero. As we have seen the parameters ( n of the gauge (12. 7h and ( 12. 8ft are polynomials in 
^--derivatives of the field ip (see appendix [B]). Therefore, for finite energy solutions we see 
that g 1 as x — >■ ±oo, and so 

n { - 1] -4 1 



4 



a, (1) 



- (oj 2 Kac. — m ) as x — >■ ±oo, (2-21) 



a| 2n+1) ^ n = l,2,. 
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We can also investigate this behaviour more explicitly by analyzing (12.121) . (I2.14p . (12. 4p and 
(I2.19p . Consequently, for finite energy solutions satisfying (I2.19p . we have that 

Of course, the theory (12.11) is invariant under space-time translations and so its energy mo- 
mentum tensor is conserved. The conserved charge is in fact a combination of the energy 
and momentum of the field configuration. In section |3] we will analyze the anomalies a^ 2n+1 ^ 
for a concrete perturbation of the sine-Gordon model, and we will show that even though the 
charges are not exactly conserved they lead to very important consequences for the dynamics 
of the soliton solutions. 

A result that we can draw for general potentials, thus, is the following. For static finite 
energy solutions the charges Q( 2n+1 ^ are obviously time independent, and as a consequence 
of (I2.22p one sees that the anomalies vanish, i.e. a^ 2n+1 - ) = 0. Under a (1 + l)-dimensional 
Lorentz transformation where x± — > x± one finds that the connection (I2.2p does not really 
transform as a vector. However, consider the internal transformation 

A^-fA^f* (2.23) 



where d is the grading operator introduced in ( 1A.3[) . Then, one notices that A^, given in (12. 2 p 



transforms as a vector under the combination of the external Lorentz transformation and the 
internal transformation (I2.23p . For the same reasons the transformed connection a M , defined 
in (12 .7p . is also a vector under the combined transformations. Consequently, the anomalies 
(3 <y2n+1 \ introduced in (I2.15p . are pseudo-scalars under the same combined transformation. 
Therefore, in any Lorentz reference frame the integrated anomalies a^ 2n+1 \ defined in (I2.18p . 
satisfy 



a 



(2n+i) _ q f or an y s ^ a ^j c or a travelling finite energy solution (2.24) 



where by a travelling solution we mean any solution that can be put at rest by a Lorentz 
boost. Even though this result may look trivial, it can perhaps shed some light on the nature 
of the anomalies a^ 2n+1 \ In fact, as we will see in our concrete example of section [3j the 
anomalies vanish in multi-soliton solutions when the solitons they describe are far apart and 
so when they are not in interaction with each other. The anomalies seem to be turned on 
only when the interaction takes place among the solitons. 

2.1 A second set of quasi conserved charges 

Note that we can also construct a second set of quasi conserved charges for the theories 
(12. ip using another zero curvature representation of their equations of motion. The new Lax 
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potentials are obtained from (12. 2 j) by interchanging x + with x_, and by reverting the grades 
of the generators. Then we introduce the Lax potentials 



A. 



A, 



(lu 2 V — m ) b-x — iuo — — F_i 

v ' dip 



1 

2 

\h - l -ud + ^F Q . 



dV 



(2.25) 



In this case using the commutation relations of appendix |X] we observe that the curvature of 
such a connection is 



d+A_ - d_A + + 



A + ,A. 



with 



X = --uid+p 



d 2 V 
1^ 



XF_ x + -u 



+ u V — m 



o2 dV 



F 



(2.26) 



(2.27) 



The construction of the corresponding charges follows the same procedure as in section [2j We 
perform the gauge transformation 



■x z-l 



A^^ a^ = gA^g - d^gg 



(2.28) 



with the group element being 

r oo 

~g = exp ^(.J_ n (2.29) 

_n=l 

and analogously to the case of section (2], we choose the C-«' s to cancel the F_„'s components 
of a + . We then have 

d + d- -d-d + + [5+, o_] =XgF_ l g- 1 (2.30) 

where we have used the equation of motion (12.11) to cancel the component of F + _ in the 
direction of Fq. The details of the calculations are given in the appendix O The transformed 
connection takes the form 



1 



a+ = - 6i + J2 s 



~(-2n-l) 



ra=0 



0-2n-l, 
oo 



~(-2n-l) k _|_ ^(-n) 

n=0 n=2 

The transformed curvature f!2.30j) leads to equations of the form 



d+ ~(-2n-l) _ d _~(-2n-l) = ~^2n-l) 



n = 0,l,2,... 



(2.31) 
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with J 5( 2n+1 ) being linear in the anomaly X, given in ( I2.27p . and the first three are given by 

= 0, 
/3 ( ~ 3) = iudiipX, 



(3 



(-5) 



°-u 2 (d + p) 2 d 2 + p + dip 



X. 



Following the same reasoning as in section [2} we find that for finite energy solutions we have 
the quasi conservation laws 



dt 







2n-l) 



dt 



-~a^-V n = i, 2 ,... 



with 



Q 



(-2n-l) 



5 (-2n-D 



(2.32) 

dxft- 2n ~ l \ (2.33) 



3 The expansion around the sine-Gordon model 

The construction of quasi conserved charges of section [2] was performed for a very general 
potential, and no estimates were done on how small the anomaly of the zero curvature con- 
dition really is. We now turn to the problem of evaluating the anomalies a^ 2n+1 \ introduced 
in (12.18p . and to discuss the usefulness of the quasi conservation laws H 2 . 2 2 1) . In order to do 
that we choose a specific potential which is a perturbation of the sine-Gordon potential and 
that preserves its main features like infinite degenerate vacua and the existence of soliton-like 
solutions. So we consider the potential given in ( 11. ip and we put n = 2 + e i.e. we take 



V(<p,e) 



(2 + e) 



tan ip 



sin cp 



\2+e 



(3.1) 



In order to analyze the role of zero curvature anomalies we shall expand the equation of 
motion ( 12. ip for the potential (13. ip . as well as the solutions, in powers of e. We then write 



ip = ip + ipi e + (f 2 e + ... 



(3.2) 



and 



dV 

d<~p 



dV 

dip 

dV 
dp 



e=0 



e=0 



+ 



+ 



d (dV s 



de \dp 



e + 



6 = 



d 2 V d 2 V dp 



dedp dp 2 de 



e + 



J £=0 
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Using the results of appendix 01 where we give the detailed calculations of such expansion, 
we have that the order zero field (fo must satisfies the sine-Gordon equation, i.e. 



d 2 ip + - sin (4<p ) = 0. 



(3.3) 



On the other hand the first order field <pi has to satisfy the equation 

d 2 ifi + cos (4</?o) <£i = sin(</? ) cos(y2o) 2 sin 2 ip In (sin 2 (<y? )) + cos 2 ((^o) • (3.4) 

We shall consider here only the anomalies for the charges constructed in section [2] (the anal- 
ysis for the charges constructed in section I2TT1 is very similar). We expand the anomaly X 
introduced in ( \2A\i as 

X = X + X 1 e + X 2 e 2 + ... (3.5) 
and we also expand the parameters 



Then we find that 



m 



IUJQ 



uq + uoi e + tu 2 e + . . . 
mo + mi e + m,2 e 2 + . . . . 

d 2 V 



|e=0 +^0 V \e=0 -m 



(3.6) 



(3.7) 



Using the results of appendix Owe find that Xq vanishes by an appropriate choice of param- 
eters, i.e. 

Xq = when uq = 4 and mo = 1. (3-8) 

With such a choice the first order contribution to X reduces to (again using the results of 
appendix [D]) 



X\ = i2d-(fo —6 sin 2 ip In (sin 2 </? ) ~ cos 2 (fo 



?7li 



+ 



y - 1 ) sin 2 tp cos 2 ip 



(3.9) 



and so we see that X\ does not depend upon tp 1 . 

Since the anomalies a <y2n+l \ introduced in (I2.18p . are linear in X and since Xq = 0, it 
follows that their zero order contribution vanishes, as it should since sine-Gordon is integrable. 
Thus we write our anomalies as 



a 



(2n+l) 



ot 



(2n+l) 
1 



e + a 



(2n+l) 2 



+ . . . 



and the first order contribution to the first two of them are (remember that a*- 1 -* = 0) 



(3) 



a 



(5) 



iu 



/oo 
dx Xi dtifo, 



(3.10) 



(3.11) 



-oo 
roc 



/oo 
dx X\ 
-oo 



- uj 2 (<9_y9 ) dlipo + ditpo 
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with X\ given in ( 13. 9p . Thus, the first order anomalies do not depend on the first order field 
ipi. The first order charges, however, do depend upon ip%. To see this we expand the charges 
as 

Q (2n + l) = g(2n+l) + g(2n+l) £ + Q (2n+1) £ 2 + ^ 

Then we find that Qg 2n+1 ^ are conserved and correspond to the charges of the sine-Gordon 
model, and involve (fo only. As an example we present the first charge at first order 

Qi 3) = f dx 8 (d-y^o) 3 (u)i d-Lpo + 4<9_y?i) + d^(p (ui d-(p + 2d_<^i) 

J — oo L 



+ 2 dt(pi d-ifQ + - sin (4 (p ) [uid_ip + 2d_(fi) - 2 d__ip d + d-ifi - i Xid^ip 

which, indeed, does depend on tpi. 

We can now evaluate the anomaly, to first order, for some physical relevant solutions of 
the theory (12.11) with the potential given by ( 13. ip . As we have stressed this earlier the first 
order anomaly depends only upon the zero order field tfo which is an exact solution of the 
sine- Gordon equation ( 13. 3p . 



3.1 Anomaly for the kink 



First we look at the case of one kink. The kink solution is given by (11.2j) with n = 2 + e and 
it is an exact solution of (12. ip for V given by (13. ip . The first order anomaly depend upon the 
kink solution of the sine-Gordon equation ( 13. 3 p which is given by 

(3.13) 



ipo = arctan (e x ) . 
Inserting this expression into ( 13.1 ip and ( 13. 91) we find that 



(3) 



a 



(5) 



dx 



sinhx 
cosh 4 x 



1 e a 
6e x ln[- 



2 cosh a; 



+ e" 



This expression can be integrated explicitly using the fact that 



d 

dx 



2sinh(x) + e^(e^-3)ln(|^) 



and so 



2 cosh'' (x) 



00 sinh x 
dx - r - 

-oo cosh X 



sinh x 
cosh 4 x 



1 e x 
6e x lnf- 



2 coshx 



+ e 



6e x In 



1 



,2 cosh a; 

Therefore, the first order anomalies vanish, i.e. = 
result shown in (12.24ft . 



+ e~ x 



(5) 



(3.14) 



= 0. (3.15) 
0, agreeing with the general 
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3.2 Anomalies for the 2-soliton solutions 



3.2.1 The soliton/anti-soliton scattering 



Let us consider a 2-soliton solution corresponding, for r\ = 1, to a soliton moving to the right 
with speed v and located at x — —L at t — 0, and an anti-soliton moving to the left with 
speed v and located at x = L at t = 0. For rj = — 1 the roles of soliton and anti-soliton 
are interchanged. The solution at order zero in the e-expansion is given by a solution of the 
sine- Gordon equation ( 13. 3ft given by 

rjv cosh?/i 



ArcTan 



sinh T\ 



with 



x 



vt- L 



n 



+ rjlnv. 



(3.16) 



(3.17) 



Putting this expression into (13.111) and (13. 9p we find that the first anomaly at first order is 



(3) 



x 



(1 - V 2 ) 3 / 2 



sinliTi coshri 



oo I 

dx — 



((3 + v 2 ) Qi + 4 v 2 ) cosh 2 yi - 2 v fii 



x 



-6 



ir cosh / v cosh \ sinh t\ 



Ai 



In 



Ai 



Ax 



— m 



2 



f 2 cosh 1/1 sinh ri 



Ax 



A: 



(3.18) 



(3.19) 



where we have introduced 

Ax = sinh 2 T\ + v 2 cosh 2 yi, Qi = sinh 2 T\ — v 2 cosh 2 yi. 

Note that given in ( I3.18p . is an odd function of T\ due to the term sinhri in front of the 
integral. All other terms involving t\ in ( I3.18P appear as cosh t\ or sinh 2 T\ , and so are even 
in T\. Consequently we see that 

[°° dtaf 1 = 0. (3.20) 
J — 00 

We point out that this result is independent of the values of uj\ and m\ which appear in the 
expression for a\ . Note that, from ( 12.221) . (13 . 101) and ( 13.121) . we have that 



dt 



(3) 



1 (3) 

"2 ai 



(3.21) 



and so 

Qf (t = 00) = gS 3) (t = -00) . (3.22) 

Thus, in the scattering of the soliton and anti-soliton the charge at first order is conserved 
asymptotically. From the physical point of view that is as effective as in the case of the 
integrable sine-Gordon theory. The solitons have to scatter preserving higher charges (at least 
in first order approximation). 
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3.2.2 The soliton/soliton scattering 



Next we consider a 2-soliton solution corresponding, for r\ — 1, to a soliton moving to the 
right with speed v and located at x = —L at t = 0, and another soliton moving to the left 
with speed v and located at x = L at t = 0. For ry = — 1 the roles of soliton and anti-soliton 
are interchanged. The solution, at order zero in the e-expansion, is again given by a solution 
of the sine-Gordon equation ( 13.31) . namely 



ArcTan 



-i] coshr 2 
v sinhy 2 



where 



x 



2/2 



VT 



+ r\ In v, 



vt- L 



T 2 



(3.23) 



(3.24) 



Following the same procedure as in the case of soliton/anti-soliton solution, by putting (I3.23f) 
into (13. lip and (13.91) we find that the first anomaly, at first order, is 



(3) 



(1-v 2 ) 



2^/2 



sinliT2 coshr 2 



An 



cosh 2 t 2 ^ / cosh 2 r 2 " 

11 V A 2 , 



w 2 sinh 2 y 2 



v 2 + 3) fi 2 -4w 2 ) sinh 2 t/ 2 + 2wfi 2 

cosh 2 t 2 w 2 sinh z y 2 



m 2 + 



^-1 
2 



A 5 



A 9 



(3.25) 



with 



A 2 = cosh 2 r 2 + v 2 sinh 2 y 2 , fi 2 = cosh 2 r 2 — v 2 sinh 2 y 2 - (3.26) 

Again, one notices that given in (13.251) is odd in r 2 . Indeed, except for the factor sinhr 2 in 
front of the integral, all other terms are even in r 2 since they involve only cosh r 2 . Consequently, 
we again have 



dt a>i 



(3) 







(3.27) 



and such result is independent of the values of ui\ and vri\. Again, using (" 13 . 2 1 f) we see that 



(t — °°) — Qf (t = —oo). So, the solitons scatter preserving higher charges asymptoti- 
cally, like in the case of soliton/anti-soliton scattering discussed above. 



,(3) 



3.3 Anomalies for breathers 

As we show in the next section the theory (12. ip with potential (13 .ip has long lived breather-like 
solutions. Hence, next we evaluate the anomaly, to first order, for such a solution. For that we 
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need the solution for the zero order field ipo which is a breather solution for the sine-Gordon 
equation (13. 3ft . i.e. 

y/l — u 2 sin(z/£) 



(po = arctan 



cosh (VT 



u 2 x 



(3.28) 



where u is the frequency of the breather (0 < u < 1), and we have chosen to express it in 
its Lorentz rest frame. Putting this configuration into (13. lip and (13. 91) we find that the first 
anomaly, to first order, is 



(3) 



4 z/ 3 (l - u 2 ) sin (2ut) h (u, t) - m\ I 2 (u, t) + 8 - l) J 3 (u, t) 



with 

h (f , t) 



dx [2 (l - u 2 ) sinh 2 (y/l - v 2 x) Q 
+ cosh 2 (Vl -u 2 x) ((l - 2 v 2 ) Q - 4 v 2 (l - z/ 
1 



(3.29) 



(3.30) 



x 



A 3 



(1 - z/ 2 ) sin 2 (yt) z 7 (1 - v 2 ) sin 2 (y t) \ ^ 2 cosh 2 (y 7 ! - v 2 ■ 



A 



A 



A 



and 



7 2 (i/,t) = y°° [2 (l - z/ 2 ) sinh 2 (Vl - ^ x) Q, 



+ cosh 2 ( y/l - u 2 x) ((l-2z/ 2 ) fi-4z/ 2 (l-u 2 



and 



I 3 (u,t) = J°° dx [2 (l - z/ 2 ) sinh 2 (Vl - v 2 x) Q 



+ cosh 2 ( y/l - v 2 x ) ( ( 1 - 2 z/ 2 ) J] - 4 z/ 2 ( 1 - v 



X 



A 3 



z^ 2 (1 — u 2 ) sin 2 (i/t) cosh 2 (y/l — v 2 x 



A 2 



(3.32) 



where we have denoted 

A = v 2 cosh 2 (y/l - v 2 xj + (l- v 2 ) sin 2 (ut) , 
Q = v 2 cosh 2 (Vl - v 2 x)-(l- v 2 ) sin 2 (ut) . 



(3.33) 

Note that the time dependence of the integrals Ij (u,t), j = 1,2,3, comes only through the 
factor sin 2 (u t) = i [1 - cos (2 ut)}. Since these integrals are multiplied by the factor sin (2 u t) 
in (13.291) . we conclude that a\ is periodic in time with period T = ^. In addition, we observe 
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that Ij (v,t) = Ij (v, — t), and so (t) = —cti (—t), due to the overall factor sin(2z/t) in 
(I3.29p . Consequently, we have that 

dt' af (t') = / dt' af ] (0 = 0, (3.34) 

J-T/2 

where we have used the fact that // +r = Jf T ^ 2 + jfr/2 + St/2: an d so the fi rs t an d third 
integrals cancel due to the fact that a\ (t) — a\ (t + T). Consequently, we find from (13.211) 
that the charge (to first order) is periodic in time 



4 Numerical support 

To check our results on the anomaly we have decided to perform various simulations of the 
Bazeia at al model - studying two kinks, a kink-antikink, and a system involving two kinks 
and an antikink. 

In all our numerical work the time evolution was simulated by the fourth order Runge - 
Kuta method. We used various lattice grids (to make sure that our results were not contam- 
inated by any numerical artefacts, the issue here was the size of the lattice and the lattice 
step). We found that to have reliable results the lattice grid (given that the kinks were of 
size ~ ±5) had to stretch to, at least, ±50. Hence most of our work was performed using 
even larger grids and the results given in this paper were obtained in simulations in which the 
lattice contained 10001 equally spaced points and stretched from -75 to 75. At the edge of 
the grid (in practice from -71 to -75 and from +71 to +75) we absorbed the kinetic energy of 
the fields. During the scattering process there was some radiation sent out towards the edges 
of the grid and it is this radiation that our procedure absorbed (so that we would not have 
any reflection of the radiation from the boundaries). Thus our procedure had the effect of 
simulating an infinite grid in which we looked only at the fields in a finite region. Thus, due 
to this absorption, the total energy seen in our simulations would decrease but this decrease 
could be associated with the system radiating some energy towards the boundaries and the 
energy seen by us corresponded to the energy of the system that we have tried to describe. 



4.1 Kink-kink interactions 



First we looked at the interaction between 2 kinks. To study this we placed two kinks at some 
distance from each other and then performed a simulation to see what happens. 
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Figure 4: Trajectories: a) n = 2, b) n = 1.9, 



When we performed this simulation with static kinks we have found that the kinks repel. 
To decide what happens during the scattering we decided to plot the positions of the kinks 
as a function of time. There are several possible definitions of 'the position of the kink' but, 
physically the most sensible one, involves looking at the energy density of each kink, with the 
position being defined at the location of the maximum of this density. This is the definition 
we have used in our analysis. 

In fig. 4 we present the trajectories of two kinks (given the definition of the position as 
mentioned above), initially at rest, as seen in simulations for n = 2 and n = 1.9. The kinks 
were initially placed at d = ±7.0 and so far away from any boundary. 

It is clear from these plots that the kinks repel. We have repeated our simulations for 
various values of n and each time the situation was the same. Looking at the plots of the 
trajectories we do not see much difference between n = 2 and n = 1.9. 

Next we sent the kinks towards each other with some velocities. In the fig. 5 we present the 
trajectories of kinks sent towards each other with velocity v = 0.5 for the cases corresponding 
to n = 2 and n = 1.9. Initially the kinks were placed at d = ±15.5 

In both cases the kinks clearly come very close towards each other but then they repel and 
move away from each other. At larger velocities they get closer before their repulsion always 
sends them back. Similar results were obtained for other values of d, n and their velocties. 

Incidentally, the sine-Gordon model (i.e. the model with n = 2) possesses a solution 
describing two moving kinks and our results (for n — 2) reproduce them very well and, 
surprise, surprise, the model does not have any static solutions involving more than one kink. 

On the other hand, the sine-Gordon moving kinks solutions are known in an explicit form, 
and because these kinks are described by explicit functions it is often said that the "kinks 
pass through each other" . This is clearly wrong when one looks at the energy density of the 



20 



20 40 60 80 



20 40 60 80 



Figure 5: Trajectories: a) n = 2, b) n = 1.9, 



moving kinks as they move towards each other. In fact, one easily observes that the kinks 
never come on top of each other (i.e. the two peaks of the energy density never form a double 
peak); in practice, the functions which describe each kink switch after the kinks' interaction. 

We have also looked at the scattering of two kinks from the point of view of the integrability 
discussed in the previous section - i.e. from the point of view of the anomalies. 

To do this we considered the scattering of two kinks for values of n close to 2. We looked 
at various positions of kinks and various velocities. All results were qualitatively similar so 
here we present our results for v = .5. The kinks were initially placed at ±15.5. We performed 
many simulations of the dynamics of such systems. In each mentioned above, the 

kinks came close to each other, reflected and then moved to the boundaries with essentially 
the original velocity. Thus the scattering was very elastic. Looking at the scattering in more 
detail it was easy to see that, strictly speaking, there was also some radiation emitted during 
the scattering and that the amount of this emitted radiation increased with the increase of 
|e| = | n — 2|; however, even for n — 1 this radiation constituted less than 2% of the total 
energy. Hence the scattering was very elastic. 

We have also looked at the values of the first anomaly and its time integrated value for 
these scatterings. 

In fig 6. we present a representative selection of our results. Fig. 6a and 6b present 
the time dependence of the anomaly and its time integrated form for n = 2, as seen in our 
simulations. Of course we know that for n = 2 the anomaly vanishes so our results provide the 
test of our numerics. We note that our values of the anomaly are very small - i. e. consistent 
with zero. Next we looked at the values of n ^ 2 for which the anomaly does not vanish. In 
fig 6c and 6d present our results for the anomaly and its integrated form for n — 1.99 and fig 
6e and 6f present similar results for n = 2.01. Fig 6g and 6h refer to the case of n = 1.98 
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while fig 6i and 6j give the results for n = 2.1 and fig 6k and 61 for n = 3. 

From these results we see very clearly that for all values of n (with the exception of n = 3) 
the integrated anomaly is approximately zero. This supports our analytical results and it 
shows that (for small e) the unintegrated total anomaly is approximately proportional to 
e = n — 2. This is supported further by the observation that the anomaly changes sign as 
e — > —e. The second order (in e) effects are comparable to those of the first order and the 
expansion in e clearly does not converge for e ~ 1. This last point is very clear from the 
case of n = 3 in which case e — 1. Of course it would be nice to understand why all the 
terms in the e power series expansion are comparable in magnitude; at this stage we have no 
understanding of this fact. 

4.2 Kink antikink scattering - quasi-breathers 

Next we have looked at the kink - antikink configurations and breathers. In the sine-Gordon 
model we do have breathers and their analytical form is well known. They are in fact bound 
states of a kink and an antikink. This is all well known; what is perhaps less known, is that 
one can generate breathers by taking a kink and an antikink and placing them not too close to 
each other and then let the configuration evolve in time. As the kink and the antikink attract 
they move towards each other, alter their shape and, at the same time, emit some radiation 
and become a breather. Interestingly, they do not annihilate but do form a breather. If we 
then absorb the energy at the boundaries the system stabilises and essentially stops emitting 
further energy as the fields have taken the shape of a breather which is a time dependent 
solution of the model. 

It is sometimes thought that the existence of breathers and of other similar configurations 
(wobbles etc) is, at least in part, associated with the integrability of the sine-Gordon model. 
Actually, as we have stressed this before, the models of Bazeia et al [lj do not appear to 
be integrable for any n other than 2; so we have decided to apply our procedure to look at 
configurations of a kink and an antikink for other values of n. However, before we discuss 
some of our results obtained in such cases let us first present them for the sine-Gordon model. 

In fig. 7 we present the time dependence of the energy of the field configuration which 
involved a kink and an antikink initially placed at d = ±3 (fig 7a) and d = ±5 (fig 7b). In 
fig 7c we present the plot of the time dependence of the position of the kink (or antikink) for 
x < 0. 

We note that in both cases (and also in other cases we have looked at, for which the 
results are not presented here) we have a very small initial drop of the total energy of the 
configuration to the energy of the resultant breather. The final energy of the breather is given 
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Figure 6: Anomalies and the corresponding integrated anomalies (from left to right 
and then down) for n = 2.0, n = 1.99, n = 2.01, n = 1.98, n = 2.10 and n = 3.00 
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Figure 7: (ab) Time dependence of total energies; (c). Position of the kink. 



by E = 2E \/l — uj 2 and its frequency u is related to the intial extend given by d. 

Next, we repeated the same procedure of generating breathers for field configurations 
corresponding to other values of n. In fig 8. we present our results for two values of n, namely 
n — 1 and n = 3.1, in which the kink and the antikink were initially placed at d — ±4.0. Our 
plots give the time dependence of the total energy of the configuration (after the absorption 
at the boundaries has eliminated the radiation reaching the boundaries). 

We note a fundamental difference; for n — 1 the energy seems to 'stabilise' around some 
finite nonzero value while for n — 3.1 it quickly goes to zero. 

We have performed many simulations (for different values of n and for different distances 
between kinks and antikinks) running them for very long times. We have found that for some 
values of n the fields annihilate very quickly; while for the others the fields evolved towards 
breather-like configurations. This was not much dependent on the distance between the initial 
kinks but depended much more on n. In fact, as the distance d increased the whole process, 
like for the sine-Gordon model, was slower, the initial radiation was smaller and the generated 
breather was larger (and so its oscillations were slower). Looking at the dependence on n it 
was clear that the closer n was to 2 the more stable the breather was (this was true from 
about ~ 0.8 to around ~ 2.8); and for some values of n (very close to 2, like 2.01 or so) the 
resultant configuration was almost indistinguishable from a breather. In fact, in all such cases, 
the energy kept decreasing but this decrease was infinitesimal. Thus we could say that we had 
a quasi-breather (i.e. a long-lived breather). As the lifetime of such a quasi-breather could 
be counted in millions of units of time, such fields, for practical (but not purely mathematical) 
reasons were not very different from a breather. 
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Figure 8: Time dependence of the total energy (a) n — 1, (b) n — 3.1. 



In fig. 9 we present some results of our simulations (for n = 2.01) which demonstrate the 
existence of our quasi-breathers. In fig 9ab we present the plots of the field configuration for 
two values of t, namely t = 355500 and t = 356300. We see that the fields look very much 
like those of the n = 2 breather. 

In fig. 10 ab and c we present the time dependence of the energy of the configuration 
on t, a detail of this dependence at large t and the time dependence of the value of field at 
x — 0. Note the extremely large values of t in the plots of the energy density. Note also the 
irregularity of the energy decrease. The energy gradually appears to decrease less and less 
and then suddenly drops and changes its slope of decrease. It then continues in the same way 
until the slope changes again etc. We do not understand these changes but, in any case, the 
total decrease of the energy is still very modest and it is clear that the quasi-breather is not 
going to "die" soon. 

The plot in fig. 10c presents the time variation of the field at x — 0. As the field is 
symmetric around x = this plot demonstrates the frequency of the oscillation. 

We have performed similar simulations starting with the initial kinks and antikinks at 
other distances from each other and for other values of n. The results were qualitatively the 
same; the further the initial structures were - the slower was the decrease in energy (i.e the 
longer the life-time of the breather). The same was true when we considered n further away 
from n — 2. This once again suggests that the models for n ^ 2 (but close to 2) are quasi- 
integrable as discussed in the previous sections. Hence we have also looked at the behaviour 
of the anomaly for our quasi-breathers. 

In fig. 11 ab we present the plots of the anomaly and the time integrated anomaly (at 
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Figure 9: Field configurations (of n = 2.01) for (a) t = 355500 and (b) t = 356300 
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Figure 10: The time dependence of the total energy (a), of the details of this 
dependence at the larger values of t (b) and of the values of the field at x = seen 
in a similation of a kink and an antikink in the n — 2.1 model (c). 
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Figure 11: The first anomaly, b) the integrated first anomaly for n = 2.01, c) and 
d) the same for n = 2.7. 



later times) of our n = 2.01 quasi-breather and in fig. 11 cd the similar plots for n = 2.7. 

We see that, when compared to the two-kink case, the time integrated total anomaly is 
nonzero at all times (as the kink and antikink are bound into a pulsating quasi-breather) but 
when one looks at the time integrated anomaly over a multiple of the period of its oscillations 
its value is very very small (i.e. close to zero). For n = 2.01 this value is so small that it is 
difficult to see that it changes at all, for n = 2.7 the value is small but nonzero and the value 
changes - thus the quantity corresponding to the anomaly is, strictly speaking, not conserved. 

4.3 Wobbles 

Finally we looked at wobbles i.e., fields involving a kink and a breather. In the sine-Gordon 
model they are again well known and, in fact, one has their analytical form. Of course, as 
before, we can generate them, numerically, from field configurations involving an antikink and 
two kinks (or vice- versa). However, as these configurations have an excess of energy, which is 
emitted when an antikink a kink form a breather, this energy can be, in part, converted into the 
motion of the remaining kink (or of the breather). Hence it is much harder, by comparison 
with pure breathers (where one can exploit the symmetry of the initial configuration), to 
generate non-moving wobbles. We have performed many simulations and the resultant fields 
sometimes were static but most of the time were moving. Clearly, the result of the simulation 
depends on the excess of energy - so further the initial structures were from each other the 
more likely there were to remain static. But this, in turn, slowed down the process of the 
generation of the breather. In addition, the futher n was from n = 2 the more radiation was 
sent out by the system and more likely it was that this radiation would set in motion the kink 
or the breather. However, for n close to 2 we did manage to obtain wobbles and in the plots 
given below we show some of our results. 
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Figure 12: a) Total energy, b) Potential energy 
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Figure 13: a) t = 0, b) * = 6400, c) t = 12800 



First we present our results for n — 2, i.e. for the sine-Gordon model. 

In fig. 12 we plot the time dependence of the total and of the potential energy seen in 
the simulation involving the kink, the antikink and the kink originally located at -15.5, and 
+ 15. In the following figure we exhibit the field configurations for three values of t, namely 
for t — 0, t — 6400 and t = 12800. We note a fast decrease of the total energy over the initial 
period and then stability. The potential energy is virtually always close to 12 and then it 
decreases to just over 4 when the breather is 'breathing', i.e. when almost all its energy is 
kinetic. Note that, for the breather, the flow of energy between the kinetic and the potential 
energies is very uneven; most of the time the breather's energy is mainly potential and the 
periods over which the kinetic energy dominates are relatively short. 

Next we present our results for the case of n ^ 2 i.e. n = 2.01. In this case the energy 
continues to decrease but this decrease is very slow. In fig 14c we present the details of 
the plot of the total energy for larger values of t. We clearly see the decrease - hence the 
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Figure 14: a) Total energy, b) Potential energy, c) Total energy for large values of 
t (n = 2.01) 



breather is slowly dying but its decay is very slow indeed. And fig 15 shows the fields at some 
representative values of time (for n = 2.01). It is clearly very difficult to see any fundamental 
difference between the wobbles in these two systems (for n = 2 and n = 2.01), although 
as fig 14b shows, the time dependences of the breather oscillations in both cases are very 
different (much shorter for the n = 2.01 wobble and slowly decreasing when compared to that 
of n = 2.00). 

As we have mentioned before, in our studies we have also seen simulations in which the 
breather and the kink which form the wobble move relative to each other. In fig 16. we show 
the field configurations for one of such cases. This case corresponds ton = 1.9. 
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Figure 15: Fields at 3 values of t; a) t = 0, b) t = 6400, c) t = 12800 
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Figure 16: Fields at 3 values of t (for n = 1.9): a) t = 500, b) t = 2000, c) t = 5000. 



5 Summary 

In this paper we have made the first steps to introduce the concept of quasi-integrability 
and discussed it on the example of the models of Bazeia et al pQ. We showed that when 
the models are close to being intergrable and so can be compared with them then one can 
introduce many quantities, which in the integrable case are conserved and which, in the non- 
integrable not conserved. One can then calculate their anomalies that are responsible 

for this nonconservation, in a power series of the difference of the quasi-integrable models from 
their integrable neighbours. In the models of Bazeia et al [1] this difference is provided by 
e = n-2. 

We can then calculate these anomalies for various field configurations. We have shown 
that the anomalies are very small for many such configurations and are only significant when 
the fields describe strongly time-dependent tightly bound objects; i.e. the field configurations 
like those of breathers or wobbles. And, conveniently, the models of Bazeia et al do possess 
such configurations. 

In fact, the models of Bazeia et al, which depend on a parameter n (which when n = 2 
reduce to the integrable sine-Gordon models) have many very similar properties and can be 
used to discuss the concept of quasi-integrability. All models (i.e. for any n) have one kink 
solutions and their scattering properties are very similar. Moreover, no other analytic solutions 
of these models (when n^2) are known. 

However, the models can be studied numerically. When we studied these models for n^2 
but close to 2 we have found that the models do, indeed, possess long-lived breather-like field 
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configurations; i.e. when we have constructed breather-like field configurations and let them 
evolve they gradually emitted some energy but this process was extremely slow; and so we 
can claim that these models (for n close to 2) possess 'very long-lived' breather-like solutions. 
Their life-time is closely related to how close n is to 2 and when n < 0.8 or n > 2.8 this decay 
was relatively fast so that the 'existence' of these states cannot be taken too seriously. 

We have also looked at wobble-like states (i.e. states involving a breather and a kink) and 
the situation was found to be similar although the range of n for which such states appeared 
to be long lived was smaller. This is partly related to our construction of such states; we 
generated them all by taking initial configurations consisting of kinks and antikinks and then 
evolving them and absorbing, at the boundaries, any energy emitted by the configuration. 
For the wobble-like states, as the configuration was less symmetric the energy was emitted 
non-symmetrically and this often lead to more perturbation of the resultant (wobble-like) field 
configuration. 

Thus, in addition to supporting our studies of quasi-integrability, our numerical results 
demonstrated also the existence of long-lived breather-like and wobble-like states. 
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A The algebra 

We consider the sl(2) algebra 

[T 3 ,T ± ] = ±T ± , [T + ,T_] = 2T 3 . (A.l) 

We take the following basis for the corresponding loop algebra 

hm+i = A m (T+ + AT_) , F 2m +i = A m (T + — A T_) , F 2m = 2X m T 3 . (A.2) 
The algebra is 

[b 2m +l j i>2n+l ] = 0, 
[ F 2m+ i , F 2n+ i } = 0, 

[F 2m , F 2n ] = 0, 

[ b 2m +l , F 2n+ i ] = —2 F2(m+ra+l) > 

[b 2m +i , F 2n ] = — 2 F 2 (, m+n ) +1 , 

[ -?2m+l j -^2n J = — 2 6 2 (m+n)+l- 

We have a grading operator 

d = To 4- 2 ), 

dA 



d = T 3 + 2A^- (A.3) 



such that 

K 6 2m+ i] = (2m+ 1) 6 2m+1 , [d,F m ]=mF m . (A.4) 



B The gauge transformation (12.71) 



The first six parameters ( n of the gauge transformation (12. 7ft . determined through (I2.10p are 
given by 

Ci = 

( 2 = 
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Ca = ii(a; 3 (^) 3 + 3a;^ 3 ) 

( 4 = -tU(p^) 2 p^uj 3 + 3p^u 
6 v 

C5 = JLi (s^ 5 ^ ' 1 )) 5 + 30a; 3 ^( ' 3 )(^( ' 1 )) 2 + 40o; 3 (^( ' 2 ))V (0ll) + 15cV°' 5) 
Ce = ±i (23(<p^)V°' 2) u 5 + 40(^°> 2 )) 3 u; 3 + 145p(°'V 'V 0,3) w 3 
+ 35(^)V(°' 4 )o; 3 + 15^ ' 6 )a;), 
where yj( > n ) = d™p. 

The first few components of transformed gauge potentials introduced in (12. 7p are given 
2 



and 



+ 6i 

+ &3 
+ &5 

= &1 
+ &3 

+ h 

1 



-^ 2 (^) 2 



+ 



(w 2 V-m) 



1 2 q2 d V 1 

-a; a (/? iuo_tpX 

A ^ dp 2 y 



dV 



— icj 3 (5_^) 3 X + — u^ditpid.tpy - iud^ipdiX + -iud 2 _pd„X 



1 



16 

dV 



dp 



- -iujcP pX + -u 2 d A p , 

+ F 2 X 
+ F 3 d-X 



+ F 4 
+ F 5 
+ F e 



-u 2 {d„p) 2 x + d 2 _x 



i 



u 2 (d^ip) 2 d_X + -u 2 d„pd 2 _pX + d 3 _X 

\u\d„p>fx + lu 2 (d_p) 2 d 2 _x + \u 2 d 2 _pd^pd_x + u 2 d 3 _pd„px + dtx 

8 2 2 



+ 
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The anomalous terms of the gauge transformation (12. 7p . i.e., those that do no vanish due 
to the anomaly X introduced in (\2A\\ . are given by 

XgF 1 g~ 1 - [o+,a_] = 



+ h 



iud 2 _tpX 
-iu 3 {d^) 2 d 2 _^X + iudtyX 



+ F 2 [-d-X] 



+ F 3 

+ F 4 

+ F 5 

+ ... 



-d 2 _x 
1 



-co 



(d.cpYd^X - u 2 d_(pdiipX - d 4 _X 



-^u 2 d-ipd 2 _ipd-X - ^(d^ipfdlX - ^oj 2 (d 2 _ip) 2 X - ^u 2 d^tpd 3 _ipX - diX 



C The gauge transformation (12.281) 



The first six parameters C_ n of the gauge transformation (12.281) and (I2.29P are given by 

1 
2 
1 



C-i = 



C-2 = -^Y 



(2,0) 



C- 



i»( W V W) ) 8 + W* 0) ). 
-z(4( V 9( 1 '°))V (2 ' 0) ^ 3 + 3^ (4 ' 0) a; 



C_ 5 = (3oj 5 (^) 5 + 30u;y 3 ' 0) (^ (1 ' 0) ) 2 + 40u; 3 (^ 2 ' )) V (1,0) + 15tV 5 ' 0) ) , 

C _ 6 = JL f (23(v9 (1 ' 0) ) V {2,0) ^ 5 + 40(^ 2 '°>) V + 145v9 (1 ' V (2 '°V (3 ' 0) ^ 3 
+ 35(^ (1 ' 0) )V (4 ' 0) ^ 3 + 15^ {6 ' 0) u;), 
where Lp( n ^ = dip. 

The first few components of transformed gauge potentials introduced in (12.281) are given 

by 

+ 2 
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and 



a_ 



+ 6-1 

+ 6-3 
+ 6_ 5 

= 6-1 

+ 6_ 3 

+ 6_ 5 
1 



16 



+ 



-Uv 

?, v 



i 2. 2 i. ~ 

-uj d + ^— + -^d +V X 



^icud+tpd+X + -iud + ipX 



+ F_ 2 

+ F_ 3 

+ F_ 4 

+ F_ 6 

+ F-6 

+ ... 



-X 

-d+X 

- l -u 2 {d+y) 2 X - d%X 



-{d+^f) d+X - -u 2 d+(pdl<fX - d 3 + X 



-lu\d + ujfx - lu 2 (d + ip) 2 d 2 + x - ^dl^d+ujd+x - u 2 d 3 + ipd+ipx - d + x 

8 2 2 



The anomalous terms of the gauge transformation (12.281) . i.e., those that do no vanish due 
to the anomaly X introduced in (12. 27ft . are given by 



XgF^g - [ 5+ , o. 



6-3 

+ 6_ 5 

+ F_ 2 

+ F_ 3 

+ F_ 4 



iud + tpX 
-iu 3 {d + <p) 2 dl<pX + iujd+tpX 



-d + X 

-d\X 
1 



-^u 2 (d+ip) 2 d+X - uj 2 d + <pd 2 + <pX - d%X 
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+ F_ 
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-\u 2 d + ipd 2 + (pd + X - u 2 (d + ^) 2 d 2 + X - \u 2 (d 2 + ipf X - ^u 2 d + pd 3 + pX - d%X 
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2 C 



D The e-expansion 



dV , 

V = V £ =o+-^ £ =o £ + ■■■ 

as 

~dV_ dVdtp 
de dip de 



V L =0 + 



e=0 



+ 



d 2 V | 2 d 2 V dip dV aV <^ / V 



e 2 + 



J e=0 



(D.l) 



Analogously, we have 

dV dV 



dip 



dip 
dV 
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+ 



d dV 



de \dtp 



e + ... 



e=0 
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d 2 V d 2 V dip 
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(D.2) 



and 



d 2 V d 2 V 



dip 2 



dtp 2 
d 2 V 



e=0 



+ 



dip 2 
d 4 V 



e=0 



+ 



+ 



d fd 2 V s 



de \ d(p 2 



e + 



£ = 



d 3 V d 3 V dip 



dedp) 2 dip 3 de 



J e=0 



d^V dip d 3 V d 2 ip d*V fdip\ 2 
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(D.3) 



Calculating we have 



V\ £=0 
dV 

dj le=0 = 4 Sin(4v?o) 



- sin 2 (2<p ) = ^[1 - cos (Aip )} 
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d 2 V 
dip 2 



=o = cos (4 <£ ) 
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and 



dV 
!k 
d 2 V 
dipde 

d 3 V 
dtp 2 ds 
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d 2 V 
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d 3 V 
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(D.4) 



e=0 



e=0 — 



e=0 — 



U=0 



-isin 2 (v? ) 



2 

-isin(2v9 ) 



2 sin 2 (v9 ) log | sin(v? ) | + cos 2 ((p ) (D.5) 
8 sin 2 p log | sin(y? ) | + cos(2</? ) + 1 

1 
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;in 4 ((/2 ) (8 log 2 | sin(c^ ) | —8 log | sin 
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(D.6) 



+ cos 
+ 8 cos 

= ^_sec 4 (v? ) 



:(4(/? )(log | sin(<y9 ) | -l) 2 - 4 cos(6c/? ) log | sin(^ ) | + cos(6y? ) + 16 
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